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ABSTRACT 

NGC 3310 is a local galaxy with an intense, ongoing starburst thought to result 
from a merger with a companion galaxy. It has several known tidal features in the 
northwest and southern regions around the main galactic disc, as well as a closed, tidal 
loop emerging from the eastern side of the disc and rejoining in the north. This loop 
appears to be distinct from the rest of the shells surrounding NGC 3310 and is the 
first of its kind to be detected in a starburst galaxy. In this work, we present U BVR 
photometry to faint surface brightness levels of this debris network, and we explore 
various strategies for modelling NGC 3310's disc and subtracting its contribution from 
each region of debris. We then compare these photometric results with the GALEV 
spectral synthesis models, and find possible material from the intruder galaxy, sug- 
gesting that the recent accretion of several small galaxies is driving the evolution of 
NGC 3310. 
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1 INTRODUCTION 

NGC 3310 is a nearby disc galaxy undergoing an intense 
starburst. This system is particularly intriguing because it 
is a field galaxy which displays signs of a recent and complex 
merging history. At a distance of about 14 Mpc (Ho = 72 km 
s -1 Mpc -1 ) and with a radial extent of 16 Kpc, including 
its most far-flung stellar debris, this relatively small galaxy 
is a system in turmoil. It has an HI mass of 2.2 x 10 9 M fl 
an estimated total mass of 2.2 x 10 10 M Q iKreeel fc Sancisil 
l200ll) and a global star formation rate of 8.5 Mq yr _1 , most 
of which is occurring in a starburst ring surrounding the nu- 
cleus. There are also many star forming regi ons throughout 
the dis c, as well as young star clusters which Ide Griis et alJ 
(2003a) find to be approximately 40 Myr old. This current 
starburst is undoubtedly largely respo nsible for the galaxy's 
blue global colour of B — V = 0.32 ide Vaucouleurs et alJ 
ll99lT) . As the light from this galaxy is dominated by young 
stars, separating the contributions of the starburst and the 
underlying disc to the overall colour of NGC 3310 is difficult. 
There is much evidence that NGC 3310 has recently un- 
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dergone at l east one significant merge r. This idea was first 
proposed bv lBalick fe Heckmanl Jl98lTl . who suggested that 
a dwarf galaxy collided with NGC 3310 and that the fa- 
mous "arrow" on the western side consisted of the remains of 
the interloper. In addition to the visible optical debris, such 
as numerous shell-like featu res and the "bow and arrow" 
J Walker fc Chincarinil Il967h ■ the radial surface brightness 
profiles of NGC 3310 we re found to follow a n 7? 1 / 4 law, which 
has been interpreted bv lSmith et alJ (I1996T) as evidence that 
the main disc has been disturbed by a recent merger. Fur- 
thermore, while the nucleus has almost s olar metallicity, the 
disc i s substantially more metal poor Ij Heckman^^Balic 
1980; iPuxlev. Hawardan. fc Mountainlll990t iPastoriza et 



1993). Because of this di screpancy, it has been suggested 
fsee fde Griis et alJ (l2003aT) for a discussion) that NGC 3310 
collided with a gas rich, metal-poor companion galaxy. An- 
other possibility, given the low metallicity of the disc, is that 
this galaxy is in the process of forming a new stellar disc. 

Work bv lKreeel fc Sancisil (1200 if) suggests that the ar- 
row feature is not a distinct kinematic element, but rather 
the beginning of an extended HI tail which is found in the 
velocity range of 804 — 953 km/s. When this tail is considered 
in conjunction with the Southern HI tail, the total mass of 
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the interloper approaches almost 10 per ce nt of the galaxy's 
total mass. In I Wehner fc Gallagher! (120051) (hereafter Paper 
I), we announced the discovery of an additional arc of stel- 
lar debris. This arc is the first closed stellar loop found in 
a starburst galaxy, and its location does not coincide with 
the gaseous, HI tidal tails. With a total stellar mass of 
4 x 10 8 A/ Q this large loop supports the idea that NGC 3310 
underwent a larger merger than previously thought. The 
possibility also remains that NGC 3310 suffered numerous 
small collisions, rather than one large merger. 

There are several ways to further explore the merging 
history of this galaxy. In Paper I we defined photometric re- 
gions, and in this work we continue our study of them. Our 
goal is to explore the origins of these debris features and to 
determine whether they are shell-like structures created by 
stars pulled from NGC 3310's ow n disc, as in the weak inter - 
action model of shell formation IIThomson fc Wright! 1990). 
or whether they represent the remains of another galactic 
companion. We then can use this information to gain in- 
sight into how the merger (s) may have progressed. In § |5J 
we present UBVR photometry for all previously defined re- 
gions of debris. We discuss the significance of this debris net- 
work's photometric colours and use several different meth- 
ods to compare them to the colours of NGC 3310's under- 
lying disc. We explore disc modelling strategies for this dis- 
turbed galaxy, and present these results in § [3] In § we 
use spectral synthesis models to further identify the star for- 
mation histories contained in the debris regions themselves 
and compare this to the current ideas of the star formation 
history of the main disc. In § we then discuss the conse- 
quences of these results for the merging history of NGC 3310. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Broad-band Imaging 

For this study we make use of data from both the WIYN 
3.5m 1 and the WIYN 0.9m telescopes. The reductions for 
the WIYN 0.9m data, taken on the nights of 2003 March 
7-11, are described in Paper I. The WIYN 3.5m data were 
taken with the MiniMosaic CCD camera on the night of 
2004 April 20, under photometric conditions. For each night, 
numerous Landolt standard fields were observed in order to 
obtain photometric solutions. While a range in airmass was 
covered for all 0.9m data, standards were only taken at the 
WIYN 3.5m at the beginning and the end of the night. Thus, 
we used the standard extinction coefficients for Kitt Peak 
and solved only for the colour terms in our data. 

The MiniMosaic data were reduced using both IRAF 2 
and MIDAS 3 . In IRAF, the MSCRED package was used to 
perform the overscan correction, which was done on a line 
by line basis for these data. It was also used to create the 

1 The WIYN Observatory is a joint facility of the University of 
Wisconsin— Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatory. 

2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

3 Munich Image Data Analysis System, provided by the Euro- 
pean Southern Observatory (ESO). 



average bias and dome flat frames and for the final bias 
subtraction and flat division. The two separate chips of the 
WIYN 3.5 MiniMosaic CCD Camera were then combined 
into a single image with the MSCIMAGE task, which remaps 
the two chips onto a single coordinate system. Once the two 
chips were merged, the multiple exposures of NGC 3310 were 
combined with imcombine and bad pixel masks were used 
to remove the gap left from the separate mosaic chips. 

2.2 Ha Images 

The Ha images of NGC 3310 were obtained on the night of 
April 7, 2003, using the WIYN 3.5m telescope under photo- 
metric conditions with seeing of ~ 1" FW H M . The images 
were processed in the manner described above before the fi- 
nal Ha image was created. The variable seeing throughout 
the night led to a small difference in the FWHMs between 
the broad and narrow-band data. The FWHM of the Ha im- 
age was approximately 1.1 pixels, or 0.16", larger than that 
of the i?-band image. To account for this, we used PSF- 
MATCH to convolve the R-band data with a Gaussian to 
broaden the FWHM profile, while still conserving the total 
flux. We then calculated the flux-ratio for stars surrounding 
the galaxy in the narrow and broadband data, scaled the 
R-band down accordingly and subtracted it from the Ha. 
This effectively removed the continuum from the Ha image 
and as a result, removed most of the stars. 

This final image can be seen in Fig. 1. The regions of 
interest defined in Paper I and the SparsePak fibre configura- 
tion have been overlaid on this Ha image. Since Ha emission 
is associated with stellar complexes having ages of <10 Myr, 
Fig. 1 demonstrates how active star formation has retreated 
to the inner parts of NGC 3310, and is not progressing at 
significant levels in most of the outer galaxy. 

2.3 Spectroscopy 

The SparsePak jBershadv et alJ 120051 120041) data for 
NGC 3310 were obtained on the WIYN 3.5m in Kitt Peak, 
Az, on December 14, 2005. Our observations were made un- 
der photometric conditions with optimal seeing (0.5"). We 
used the echelle 316@63.4 grating at order 8, centred at 
A = 6655A, resulting in a resolution of 11,400 or 0.21 A 
per pixel. This corresponds to a resolution of 20 km s" 1 
over a two pixel resolving element. The fibre array was cen- 
tred on the galaxy's bulge, and the fibre positions are shown 
in Fig. 1. We obtained 2 exposures of 1800s each. We also 
took calibration arcs using the CuAr lamp throughout the 
night. 

The basic reductions, such as the bias and overscan cor- 
rections, for these data were done using IRAF with the CC- 
DRED and HYDRA packages. The combined images and 
calibration frames were created using the IMCOMBINE pro- 
cedure and subsequent was completed using the DOHYDRA 
task. Using this task, we flat-fielded the spectra and per- 
formed the wavelength calibration with the assistance of 
the NOAO calibration spectra database 4 . The spectra were 
linearised to facilitate line extraction and further analysis. 
Also, we used the SparsePak sky fibre data to remove the 

4 http:/ /www. noao.edu/kpno/specatlas/index. html 
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Figure 1. Continuum subtracted Ha image of NGC 3310 taken with the WIYN 3.5m. Regions of debris defined in Paper 
I have been overlaid in Figure la (left), as well as the SparsePak array of fibres. From this image, we can determine that 
there is little to no star formation in all regions of debris, except for those which define the "arrow" (regions K & L). Figure 
lb (right) shows more of the structure in the Ho emission, and the well known Jumbo HII region and the circumnuclear 
starburst ring are labelled. 



sky lines from the final spectra. Of the seven sky fibres, four 
were found to contain target Ha lines and were removed. 
The remaining three fibres were averaged and used for sky 
subtraction. 



3 PHOTOMETRY 
3.1 Integrated Colours 

The colours for each region are listed in Table 1, and the 
UBVR surface photometry for the network of debris in 
NGC 3310 is shown in Table 2. These calibrated magnitudes 
and their errors were calculated using the POLYPHOT 
package in IRAF and are described in Paper I. The errors 
shown include only random sources of error, and thus rep- 
resent lower limits. However, for low surface brightness pho- 
tometry, the main source of systematic error is fiat field er- 
ror. With the S2KB chip and the well-baffled WIYN 0.9m 
telescope, the signal is flat across the chip to within 0.7 per- 
cent. Thus, the variation across each region shown in Fig. 
1 is negligible for the 0.9m data. However, for the [/-band 
data, taken on the WIYN 3.5m, scattered light becomes a 
more serious issue and the field is less flat by an order of 
magnitude. As a result, the errors in U — B may be as high 
as ±0.1 mag. Another source of systematic error is that the 
natural variation in the colour within each region. However, 
average colours are still useful for studying trends as well 
as the aggregate properties of the stellar populations within 
each region. 

While there is quite a range in B — V colour across the 
system, in general, most of the defined regions have bluer 
colours t han a re typical of non-starbursting disc galaxies 
ide Jondll996t) . Ideally, we would like to compare the debris 
colours with those of NGC 3310's underlying disc. The dif- 
ficulty with this comes from the fact that when we observe 
NGC 3310, we are seeing light from both the underlying disc 
and the overlying starburst, combined. 



While modelling these two components separately is 
tricky to do, one can use several different approaches 
to estimate the colour of the underlying disc. If we as- 
sume that before the collision and its subsequent starburst, 
NGC 3310 was a typical d isc galaxy, we can use its RC3 
Ide Vaucouleurs et alJll99Jl classification (T = 4.0 ±0.3) to 
compare its colour with the av erage col o ur for similar galax- 
ies. According to the work of Ide Jonel jl996t) galaxies of a 
similar classification type to NGC 3310 have a B — V = 0.74. 
Each tidal feature in NGC 3310 is significantly bluer than 
this. Since the weak interaction model suggests that shells 
are formed from a galaxy's own disc stars, we expect that 
shells formed in this way would have stellar populations sim- 
ilar to the galaxy's main disc. Thus, their blue colour sug- 
gests it is unlikely that NGC 3310's debris was formed by 
a weak interaction. However, this assumption breaks down 
if the merger(s) that NGC 3310 experienced caused it to 
change its structural class or if NGC 3310's underlying disc, 
due to its starburst and possible recent interactions, is sim- 
ply not representative of an ordinary spiral galaxy, and thus 
its RC3 classification isn't as well correlated to its current 
colour as a normal galaxy would be. As discussed in further 
sections, this latter argument appears to be the case. 



3.2 Disc Models 

We turn to disc modelling to further explore the properties 
of the underlying disc. For this task, we used several dif- 
ferent MIDAS scripts to calculate the surface brightness as 
a function of radius. For this we used the FIT/ELL3 rou- 
ti ne in MIDAS, which is base d on the iterative algorithm 
of Bender & Mollenhoff (1987). Thus, the resulting surface 
brightnesses represent the two-dimensional flux distribution 
based on a \ 2 convergence fit. As such, these profiles include 
light from the underlying disc, the starburst, and overlap- 
ping stellar debris. Ideally we would like to separate these 
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Table 1. Colour Indices of the Tidal Debris of NGC 3310 



Region 




U 






B 






V 






R 






U-B 




B — V 




V — 


R 


A 


17 


,43±0. 


006 


17 


,49±0 


.02 


17. 


08±0. 


.003 


16. 


73±0. 


.004 


-0 


,06±0. 


02 





.41±0.02 


0. 


36±0 


.007 


B 


16 


.23±0. 


.003 


16 


.41±0 


.01 


15. 


93±0. 


.002 


15. 


53±0. 


.002 


-0 


.18±0. 


01 





,48±0.01 


0. 


40±0 


.004 


C 


16 


.01±0. 


002 


16 


,20±0 


.01 


15. 


77±0. 


.002 


15. 


39±0. 


.002 


-0 


.18±0. 


01 





42±0.01 


0. 


39±0 


.004 


D 


17 


.21±0. 


006 


17 


.55±0 


.02 


17 


12±0. 


.004 


16 


73±0. 


.005 


-0 


.34±0. 


03 





43±0.03 


0. 


39±0 


.009 


E 


17 


.08±0. 


006 


17 


.41±0 


.03 


16. 


89±0. 


.004 


16. 


49±0. 


.005 


-0 


.33±0. 


03 





52±0.03 


0. 


41±0 


.009 


F 


17 


.35±0. 


.007 


17 


.64±0 


.03 


17. 


,01±0. 


.004 


16. 


59±0. 


.005 


-0 


.29±0. 


03 





63±0.03 


0. 


42±0 


.009 


G 


17 


.46±0. 


009 


18 


16±0 


.05 


17. 


63±0. 


.013 


17 


24±0. 


.010 


-0 


.70±0. 


06 





,53±0.06 


0. 


39±0 


.017 


H 


16 


.54±0. 


.004 


16 


.74±0 


.02 


16 


31±0. 


.003 


15. 


93±0. 


.003 


-0 


.20±0. 


02 





43±0.02 


0. 


38±0 


.006 


I 


17 


98±0. 


.010 


19 


.26±0 


.10 


18. 


92±0. 


.025 


18. 


65±0. 


.023 


-1 


.28±0. 


10 





34±0.10 


0. 


27±0 


.038 


J 


16 


.22±0. 


.004 


17 


.83±0 


.05 


17. 


22±0. 


.008 


16. 


72±0. 


.007 


-1 


.61±0. 


05 





61±0.05 


0. 


50±0 


.013 


K 


16 


89±0. 


003 


17 


,48±0 


.02 


17. 


27±0. 


.003 


16 


91±0. 


.004 


-0 


.59±0. 


02 





22±0.02 


0. 


36±0 


.006 


L 


17 


.17±0. 


004 


18 


.05±0 


.03 


17. 


79±0. 


.007 


17 


39±0. 


.006 


-0 


.88±0. 


03 





,26±0.03 


0. 


41±0 


Oil 


M 


18 


00±0. 


007 


18 


18±0 


.03 


17. 


.77±0. 


.004 


17. 


.44±0. 


.006 


-0 


,19±0. 


03 





,41±0.03 


0. 


34±0 


.010 



Table 2. Photometric and Physical Properties of the Tidal Debris of NGC 3310 



Region 


W 


MB 


(S/N) pix a 


(S/N) ap b 




(S/N) plx 


(S/N) ap 


l-'R 


(S/N) pix 


(S/N) ap 


A 


24.10 


24.16 


0.28 


10.09 


23. 


.75 


1.52 


54.64 


23.39 


2.13 


76.47 


13 


23.94 


24.12 


0.29 


16.90 


23. 


.64 


1.68 


97.41 


23.24 


2.44 


142.01 


C 


23.80 


23.99 


0.33 


19.80 


23. 


.56 


1.80 


108.52 


23.18 


2.59 


155.92 


D 


24.17 


24.51 


0.20 


8.37 


24. 


.09 


1.12 


45.95 


23.69 


1.62 


66.63 


E 


24.38 


24.71 


0.17 


8.15 


24. 


.19 


1.01 


48.61 


23.79 


1.48 


71.28 


F 


24.40 


24.70 


0.17 


7.41 


24. 


.07 


1.13 


48.71 


23.65 


1.69 


72.38 


G 


24.98 


25.67 


0.07 


3.71 


25 


.15 


0.42 


22.37 


24.76 


0.61 


32.32 


H 


24.32 


24.52 


0.20 


12.12 


24. 


.09 


1.11 


66.76 


23.71 


1.59 


95.50 


1 


24.55 


25.83 


0.06 


2.06 


25. 


.50 


0.31 


10.51 


25.22 


0.40 


13.65 


J 


24.11 


25.72 


0.07 


4.26 


25 


.11 


0.43 


27.41 


24.61 


0.69 


43.74 


K 


22.79 


23.38 


0.57 


14.37 


23. 


.16 


2.60 


65.42 


22.80 


3.65 


91.84 


L 


23.33 


24.21 


0.27 


7.59 


23. 


.95 


1.26 


35.88 


23.54 


1.85 


52.57 


M 


23.88 


24.06 


0.31 


7.66 


23. 


.65 


1.67 


41.61 


23.31 


2.29 


57.12 


N 






0.00 


0.00 


25. 


.55 


0.29 


13.34 


25.44 


0.33 


14.95 



Signal-to-noise per pixel on the final calibrated image. 
Total signal-to-noisc for each aperture. 



components in order to compare the colours of the stellar 
debris with that of the underlying disc of NGC 3310. 

Once we obtained our best 2D flux distribution for 
NGC 3310, we subtracted it from our original image in order 
to examine the structure of the debris. This image is shown 
in Fig. 2a. In this image, the debris network and the clumpy 
star-forming regions are clearly visible. However, when we 
examine the resulting 2D distribution, we find that the pres- 
ence of debris distorts the fit on the western side (Fig. 2c) 
In order to correct for this, we wrote a new MIDAS script 
to use a slice of the eastern portion of the galaxy, shown 
as a black rectangle in Fig. 2b, which we then azimuthally 
expanded to represent the entire surface of the disc. This 
is only a partial solution in that it still contains light from 
some star- forming regions, but it allows us to sample the 
less optically disturbed side of the galaxy and to choose a 
slice containing no bright tidal debris. We then subtracted 
this fit from our image, shown in Fig. 2b. The actual 2D 
fit from this circular model is shown in Fig. 2d. In compar- 
ing Figs. 2a and 2b, one can clearly see the asymmetries 
induced by the western regions of tidal debris in the outer 



disc in 2a. These features disappear when using the eastern 
slice to model the disc. 

Next, we applied our photometric solutions to calibrate 
the model frames in each band, and subtracted them to cre- 
ate a B — V frame for the disc model. We then used the 
D25 = 3'.1 value defined in RC3 as our diameter and per- 
formed aperture photometry to find a total B — V of 0.30. 
This agrees with previous measurements of the global B — V 
for this system and provides a useful consistency check on 
our 2D flux distributions. While these models do not sepa- 
rate the disc and starburst contributions to the total emis- 
sion, they are effective at subtracting the outer regions of the 
disc, since at higher radii the contribution from the starburst 
becomes insignificant. The residuals from this model, which 
highlight clumpy galactic regions and regions of debris, can 
be seen in Fig. 2b. 

In order to carry out detailed surface brightness pho- 
tometry of NGC 3 310 we then derived th e BVR profiles 
using method iv of iPapaderos et"al] |2P02). These results 
are shown in Figs. 3 and 4. In Fig. 3a, the surface bright- 
ness is shown as a function of radius plotted linearly. The 
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(C) 



(d) 



o 



Figure 2. Left: (a) Residuals obtained by subtracting the 2D flux distribution (c) calculated from the FIT/ELL3 script in 
MIDAS from the B-band image. Right: (b) Residuals obtained by subtracting the circularly symmetric 2D light distribution 
(d) from the B-band image. The black rectangle in (b) indicates the slice of the Eastern side of the galaxy used to create 
the circular model. 



solid line is a fit to the B-band surface brightness profile 
(SBP) for R 60", and yields an approximation for the 
disc and debris at these radii. This fit implies a central sur- 
face brightness of /^o=21.77 ± 0.06 mag arcsec -2 , compara- 
ble t o the Freeman value of 21.65 mag arcsec -2 iFreemanl 
1970) . This value is somewhat lowe r than that found by 
B agge tt. Baeeett fc Anderson! jl99ST) (jity = 23.98) but our 
models include contributions from both the bulge and the 
disc, and can therefore be expected to be brighter. We also 
find an exponential scale length of 2.08±0.04 Kpc, or 31", 
and an absolute total B magnitude of -18.4 mag (for D = 14 
Mpc) for the disc. This sc ale length is somewhat longer than 
that of lDuric et ail Jl99rj) (I = 7.3") for an off-Ha waveband 
and significantly smaller than the I = 123" measured by 
iBaggett. Baggett fe Anderson! lll998l) in the V-band. How- 
ever, the SBPs we present in Figs. 3 and 4 extend at least 
two magnitudes deeper than those previously measured. 

Fig. 3b shows the B — V and B — R radial colour pro- 



files for NGC 3310, calculated by subtracting the respec- 
tive SBPs. The colours for the outer disc region (defined as 
R ^ 60") are B — V = 0.42 ± 0.03 and B — R = 0.86 ± 0.07 
magnitudes. 

It is also interesting the shape of the surface brightness 
profiles for each band. In Fig. 4 we plotted the SBPs against 
R 1 ' 4 in order to check to how well the SBP can be approxi- 
mated by the de Vaucouleurs law. From this plot we see that 
the de Vaucouleurs law (fit by r 1 ^ 71 where r/ = 4) describes 
parts of the SBP (between R 1/4 = 1.7 and R 1/4 = 2.7) 
very well, but not the SBP as a whole. A Sersic fit yields 
an 77 = 2.8 and a slightly better x 2 than a model with a 
fixed r\ — 4. However, all Sersic distributions leave system- 
atic residuals when fitted to the total SBP, indicating that 
a single fitting law cannot perfectly fit the SBP. In fact, an 
exponential disc appears to emerge at higher radii. 

It is clear from Figs. 3 and 4 that the light from the 
starburst makes a large contribution to the overall light 
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Figure 3. (a) Surface brightness profiles in BVR for NGC 3310. 
The grey line represents a fit to the B-band profile at radii ^ 
60". The ij-band is represented by red triangles, V by green 
open circles and B by blue closed circles. When error bars are 
not shown, the errors do not exceed the size of the point, (b) 
Radial colour profiles for NGC 3310. The B — V profile is shown 
in blue, closed circles, and the B — R profile is shown with red, 
open circles. 
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Figure 4. Surface brightness profiles in BVR, plotted versus 
R 1 / 4 . The lower grey curve is the fit from Fig. 3, while the up- 
per grey line shows a Sersic fit to the B-band profile. This fit 
yields 1/r] = 1/2.8, rather than the l/rj = 1/4 expected for a de 
Vaucoulcurs profile. 



(1996), but is rather more typical of an Sd galaxy. It is, 
however, in excellent agreement with our average colour for 
the outer radii of the surface brightness profiles. This sug- 
gests that B — V ~ 0.4 well represents the outer regions of 
the disc of NGC 3310. While debris within D 25 still adds to 
the total colour of the disc, its contribution is small since 
the debris is expected to be fainter than the galaxy itself. 
Furthermore, the photometry of Table 1 suggests that the 
debris in this system would bias the galaxy toward a redder 
B — V. Given the blue nature of the B — V resulting from 
the pixel masks, it is unlikely that overlying tidal debris is 
a significant contributor to the global galactic colour. 



from NGC 3310. In order to better understand the nature 
of NGC 3310's underlying disc, we must find some way to 
eliminate the light contribution to the bright starbursting 
regions. To accomplish this, we have created pixel masks for 
each of the B, V and i?-band images in order to remove the 
contribution of obvious star-forming regions across the disc. 
The same mask was used for all images, and was created by 
eliminating all counts greater than 5<r above the sky in the 
B-band image, which most closely represents the younger 
star-forming regions. While masking reduces the total num- 
ber of pixels available to the models, we have sufficient pix- 
els in regions outside the main starburst ring to measure 
colours with a reasonable signal-to-noise ratio; we have ef- 
fectively masked out the entire central region of NGC 3310, 
leaving only the outermost disc. Furthermore, most of the 
outer debris shown in Fig. 1 lies outside of the D25 radius 
of 92", and will not affect the final colour. The B-band im- 
age and its pixel mask are shown in Figs. 5a and 5b. When 
calculated using only pixels within D25 and not eliminated 
in the masks, we find a global disc colour of B — V = 0.44, 
which is still bluer than the expected value from Ide Jon j 



3.3 SparsePak Preliminary Results 

For each fibre in which Ha was detected, we fit the line 
with a Gaussian profile to find the peak of the line and 
measure its velocity. Fibres that contained lines with S/N > 
3 are shown in Fig. 6. All but two lines were well fit by a 
single Gaussian profile. For the two remaining fibres, which 
exhibited a double-peaked structure at the very top of the 
line, the profile was fit as a single line to obtain an average 
velocity. In both cases, each actual peak's velocity deviates 
from the average fit by less than 3 per cent. 

Since the SparsePak data traces the warm ionised gas 
component in NGC 3310, these data suggest that NGC 3310 
has maintained its disc structure and rotation in the star- 
forming component of its disc despite recent substantial 
mergers with this system. While the old stars may have 
experienced a substantial disturbance, the young stars, as 
represented by the Ha emission shown in Fig. 6, appear 
to exist in a rotating system. The velocities shown in 
Fig. 6 are not corrected fo r inclin ation. From HI data, 
iMulder. van Driel. fc Braind <1995ft find i = 52°. Earlier 
work focusing on Ha emission finds a slightly lower inch- 
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Figure 5. (a) _B-Band image of NGC 3310 taken with the WIYN 0.9m telescope on Kitt Peak, (b) Same B-Band image of 
NGC 3310 with all pixels great than 5<r above the sky background removed. This mask was used in all bands to estimate the 
colour of the outer disc. 



nation of i — 32deg. While our observed Av = 100 km s^ 1 
across NGC 3310's disc is at first glance somewhat low for a 
typical disc galaxy, when corrected for inclination (adopting 
the HI value) this An becomes 124 km s _1 , yielding a rota- 
tional velocity of 62 km s _1 . While the HI exhibits a higher 
total velocity, the HI data extend well beyond the range 
of our SparsePak data. At a r adius of 40", or 2.7 Kpc, our 
value is consistent with that of lMulder. van Driel. fc Brainel 
( 1995). It is also interesting to note that our SparsePak data 
show an initial increase in the velocity of the Ha emitting 
gas, and then drops at higher radii to values consistent with 
those f rom HI. This trend is al s o not ed in Ivan der Krufd 
Jl99ll) : iGalletta fc Recillas-Cru3 il982Tl The fact that the 
warm ionised gas is fairly smoothly rotating may provide 
some support for the idea that the disc in this system may 
be in the process of reforming in the m anner described 
bv lFirmaivi J _Ije£nanciez 1 j£_G^ll^gl2g rl il996f) after a gas-rich 
merger dSpringel fc Hernquisdl2005ir Or, rather NGC 3310 
may have maintained its disc structure despite multiple col- 
lisions. If this is indeed the case, galactic discs may be sur- 
prisingly robust. 

3.4 Discussion 

The B — V calculated for the outer disc regions is redder 
than the global B — V given in the RC3 which includes the 
bright star-forming regions. However, B — V for the debris 
regions ranges from 0.63 in region F to the South to 0.22 for 
the bright star- forming region, K, in the arrow. That sev- 
eral regions are still significantly redder than the underlying 
disc stars suggests that they did not originate in the disc of 
NGC 3310. One possibility is that they consist of the broken 
remains of an incoming, disrupted companion galaxy - as has 
been suggested for a similar debris loop around NGC 5907 
JZheng et all Il999t IShang et alJ Il99l : ISackett et all Il994t 
iMorrison et alJl994tT The bluer regions may reflect a subse- 
quent (and currently fading) star formation episode among 
the debris. 
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Figure 6. Radial velocity plot of NGC 3310. Velocities were 
measured from the Ha line in each fibre of the SparsePak array, 
and only those lines with S/N > 3 are shown. 



We also note the presence of an interesting colour gra- 
dient in NGC 3310's debris. In Figs. 7a and 7b, we show 
colour plots of NGC 3310 in B — V for two different colour 
ranges. In these images, a colour gradient is apparent from 
regions F to E to D to C, changing from B — V = 0.63 in re- 
gion F to B — V = 0.42 in region C. This could be caused by 
a propagating wave associated with an incoming companion 
inducing star formation along the outer debris in the arc 
defined by F-C. After a burst of star formation, the U — B 
colour will fade faster than B — V, which is less sensitive 
than the (7-band to the presence of the most massive stars. 
Thus, the increasingly blue B — V, yet fairly constant U — B 
from F to D, may represent a fading burst passing from F to 
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D. In this case, the arrow structure may represent the front 
edge of a propagating star-forming wave, leaving regions of 
recently compressed gas and new star formation in its wake. 
Likewise, the star formation rate in these red regions (F, J) 
might have been truncated some time ago if they consist of 
material expelled as tidal structures from within NGC 3310. 
In models of stellar populations that have had their star for- 
mation completely truncated, it takes about 2 x 10 s years 
to drive B — V from 0.42 to 0.62 after completely stopping 
star formation. This may also explain the colour gradient in 
regions FEDC. 

It remains a possibility that the masks were simply in- 
sufficient to remove all the light from the starburst. In this 
case, the pre-burst disc stars may be redder, and the shells 
could have originated in the disc. We further explore the ori- 
gins of these regions by comparison with spectral synthesis 
models in the following section. 



4 COMPARISON WITH GALEV MODELS 
4.1 Model Assembly 

The GALE V ste llar population synthesis models 

iSchulz et alJ l2002t lAnders fc Fritze-v. Alvensleben! [2003) 
involve two different components - one of which repre- 
sents the underlying, pre-burst galaxy, and another which 
represents the starburst. We can then compute colours 
for the combined system of disc and starburst. GALEV 
provides two options for the underlying disc component. 
The first is a galactic disc that has experienced constant 
star formation (for which we chose 2 Mq yr^ 1 ) throughout 
its lifetime. This model, when combined with a starburst 
model, defines the blue limit for a galaxy. The other option 
is to represent the underlying galaxy with an exponentially 
declining star formation rate. This model begins with a 
high rate of star formation that declines on a timescale 
of 10 9 years. In this model, the galaxy is dominated by 
passive evolution and represents the red limit for a galactic 
disc. Each of these models have been calculated with five 
different metallicities. For comparison with our work, we 
have chosen an intermediate metallicity value of Z = 0.008, 
since while NGC 3310 is thought to ha ve greater than 
solar m etallicity in the nucleus, results from lPastoriza et alJ 
dl993h suggests that material outside the starburst ring is 
subsolar in metallicity. 

The output for each model is a set of luminosities in 
different filters for a given stellar mass for each time step 
from t = to beyond a Hubble time. We combined the two 
models using: 

L\ t tot = L\,ug + bLx,SB (1) 

where L\ )Ug is the luminosity of the underlying galaxy 
component, per solar mass at a given wavelength, L\ : sb is 
the luminosity per solar mass of the overlying burst, and 
b is the strength of the burst by mass fraction. For L\ tUg 
we chose the luminosity at 12 Gyr to represent an existing 
disc galaxy. The luminosities (and colours) of the underly- 
ing galaxy are actually quite insensitive to age after approx- 
imately 5-6 Gyr, and thus colour is mostly dependent on 
which model is chosen for the underlying galaxy. 

The total luminosities were found for each band and 



were converted to magnitudes and then colours. These 
colours form the evolutionary tracks seen in Fig. 8. 



4.2 Model Results 

In order to explore the stellar populations resident in each 
debris region of NGC 3310, we compared our photometry 
of the regions displayed in Fig. 1 with various models from 
the GALEV spectral synthesis models. In Fig. 8, we show 
our photometric results plotted over several model curves 
created using the GALEV models. 

Each curve in Figs. 8a and 8c represents an underlying 
galaxy with exponentially declining star formation that has 
a starburst overlaid. The burst is represented as an instan- 
taneous event that takes the underlying, ancient galaxy (t = 
12 Gyr) from the reddest point in U — B and B — V space, 
to the bluest. The colour then slowly grows redder as the 
stars which were formed in the burst then age. Each curve 
represents a different strength of burst, ranging from 0.8 to 
2 per cent of the local stellar population by mass, and whose 
light then adds to that of the ancient stellar population. 

While the error bars are most likely higher in the U- 
band due to greater systematic errors, there are still trends 
evident in Fig. 8a. Regions A, B, C, H, K and M are all con- 
sistent with burst strength of 0.014, i.e. a starburst involving 
1.4 per cent of the local stellar population's mass. They are 
also consistent with a post-burst age of 50 ± 20 Myr. While 
region D is slightly bluer in U — B, it is also consistent with 
the rest in age, although it requires a slightly lower burst 
strength, 6 = 0.007 i.e. with only 0.7 per cent of that re- 
gion's total mass. It is possible that region D has simply 
experienced less star formation than the aforementioned re- 
gions. These results are consiste nt with the burst age found 
from NGC 3310 star clusters bv lde Griis et alJ <2003bft . 

While region L is lower than all the models, it is still 
consistent with the rest of the points when the errors of the 
models themselves are considered (±0.1 magnitudes). 

Figs. 8b and 8d show the results for an underlying 
galaxy with a constant star formation rate plus a starburst. 
Note that the initial B — V — 0.44 for the constant star- 
forming model at t = 12 Gyr. This is the same value of B — V 
that we measure when we removed the starburst component 
from the disc light using pixel masks. If these masks success- 
fully removed the starburst contribution to the total light, 
then NGC 3310 should be well represented by a constantly 
star-forming disc before the most recent burst occurred. 

In comparing our photometry to the models in Fig. 8b, 
we find that while a few of the bluer regions, such as K, 
A, M, C, and H overlap with these evolutionary tracks, the 
rest of the debris network is inconsistent with these models, 
even when the model errors of ±0.1 are taken into account. 
While it is possible these regions have a completely distinct 
evolutionary history from the rest of the debris, it is more 
likely that their evolution is connected with the rest of the 
debris system, and that their bluer colours indicate higher 
burst strengths in these regions, as is suggested by Fig. 8a. 
We therefore conclude that the exponentially decaying star 
formation rate plus burst model gives the better description 
for these data. 
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Figure 7. B-V colour map of NGC 3310. B- and V-Band images were smoothed using a Gaussian profile to match their PSFs, calibrated 
photometrically and then subtracted. The image on the right shows a smaller range in colour than the B — V range shown in the left. The 
locations of regions F, E, D, and C are shown in each image. 
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Figure 8. Left: (a) Lines show the models of exponential star formation with various burst models overlaid. The burst models range 
in strength from 0.08 to 2 per cent of the local stellar population by mass. Each of these bursts is represented by a curve, while lines of 
constant time connect the burst models. The figure below (c) shows an enlarged subset of the top figure. Right: (b) Same as the left, 
except underlying galaxy is represented by a constantly star-forming disc galaxy. This disc has a colour of B — V = 0.44, chosen to be 
similar to the average colour of the outer regions of NGC 3310. The bottom image (d) shows an enlarged subset of the figure directly 
above. 
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5 SUMMARY 

We have explored the photometric properties and the origins 
of NGC 3310's tidal debris. This galaxy is currently under- 
going a starburst induced by a complex merger. We find a 
global B-V = 0.30 for NG C 3310 which is consistent with 
Ide Vaucouleurs et alJ l)l99ll) . In an attempt to separate the 
contributions of the starburst and the underlying disc in the 
NGC 3310 system, we created a mask in the B-band to elim- 
inate all luminous blue regions in the disc. This resulted in 
a global B — V = 0.44 for areas with little active star forma- 
tion. This is also consistent with the average colour (B — V 
= 0.42) calculated from the surface brightness profiles at the 
outer radii of NGC 3310. However, the colours of the shell 
and debris structures range in B — V from as red as 0.63 to 
as blue as 0.22, as measured for region K in the star-forming 
arrow. This suggests that the debris did not originate in the 
disc. 

Another possibility is we were unable to completely sep- 
arate the contribution of the bright starburst and the un- 
derlying disc, and that NGC 3310's underlying disc is truly 
redder than B-V = 0.44. Since NGC 3310 is UV bright any 
little remaining contamination from the burst could skew the 
results to a bluer B — V. In order to explore this, we assumed 
a red, underlying disc comprised of an aged stellar popu- 
lation experiencing exponentially declining star formation 
and overlaid starburst models of various strengths to create 
composite colours. We also created another set of composite 
colours using a bluer underlying disc model with a modest 
yet constant star formation. The model for the underlying 
galaxy in this latter case also matched our calculated under- 
lying disc colours of B — V = 0.44. 

The comparison of colours of debris structures with pre- 
dictions of GALEV stellar population synthesis models also 
suggests that these structures are not simply remains of 
NGC 3310's pre-merger stellar disc. They are more likely 
the remains of a redder, disrupted companion galaxy that 
was destroyed in merging with NGC 3310. In comparing the 
B — V colours of regions in the debris with those of models, 
we find that they are consistent with an ancient stellar pop- 
ulation that has experienced a weak burst of star formation 
in the past 50 Myr. In general, the distinctive outer features 
of NGC 3310 have colours that imply significantly different 
evolutionary histories than those of its main disc, and are 
not consistent with the bluer constantly star-forming disc + 
starburst models. 

The extensive debris network - specifically the large stel- 
lar loop, the two extended HI tails and the stellar shells - as 
well as the blue underlying disc colours, the high gas content 
and the starburst may be indicative that we are witnessing 
the result of a major merger or multiple, smaller, simultane- 
ous mergers and that the disc of NGC 33 10 is being reformed 
as suggested bv lKregel fe Sands! <l200ll) . 

Our deep surface photometry indicates that the SBP 
of NGC 3310 is best approximated by a Sersic model with 
an exponent 1/r) = 1/2.8, slightly larger that the value 
I/77 = 1/4 that is expected for the de Vaucouleurs law. How- 
ever, systematic residuals between the Sersic fit and the ob- 
served SBPs suggest that a single fitting law cannot provide 
a perfect fit to the SBPs. This fact, in connection with the 
exponential slope observed in all bands for R > 60" suggest 
that a disc+starburst model gives a better approximation to 



NGC3310's total light and supports the idea that the disc 
has survived the recent mergers experienced by NGC 3310. 
T his does not entirely rule ou t a major merger, however, 
as ISpringel fc H crnquist (2005) suggest that if disc galax- 
ies are sufficiently gas rich, then discs can survive equal 
mass mergers, implying that discs may be surprisingly ro- 
bust. However, since there is some evidence of older stellar 
populations in this system, it is perhaps more likely that 
we are seeing the growth of NGC 3310's disc through the 
accretion of multiple, small companions onto a preexisting 
disc. 

Most of the debris is consistent with a small burst of 
star formation in the past 50 Myr, (which is consistent with 
the ages of the young stellar cluster s in the central star- 
bursting regi o n ( j de Griis et al . 2003a)), and recent work by 
iKnapp et all i2006l) suggests that a secondary star forma- 
tion event, as indicated by the presence of a young stellar 
cluster, may have occurred in this system approximately 750 
Myr ago. This discrepancy may result from either unique, 
sequential mergers with NGC 3310, or the incoming com- 
panion and its resulting debris completing multiple orbits 
before finally being destroyed. Either case may help account 
for the observed range of colours in the outer debris network. 

Future 2D integrated field unit spectroscopy with im- 
proved velocity resolution would be useful for further exam- 
ining the stability of the disc and exploring the likelihood 
that we are witnessing a reforming disc. It may also be use- 
ful to obtain kinematic information on the shells, although 
their low optical surface brightness and low HI content will 
make it difficult to obtain reasonable signal-to-noise ratios 
in these regions. In lieu of such data, dynamical modelling 
of this system would provide further insight into the origins 
of the surrounding stellar debris. 
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